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Abstract 



o 



^vq In a companion paper, a reduced model for propagation of acoustic waves in a cloud of inertial cavitation bubbles was 
■proposed. The wave attenuation was calculated directly from the energy dissipated by a single bubble, the latter being 
^ estimated directly from the fully nonlinear radial dynamics. The use of this model in a mono-dimensional configuration 
I 'has shown that the attenuation near the vibrating emitter was much higher than predictions obtained from linear theory, 
and that this strong attenuation creates a large traveling wave contribution, even for closed domain where standing waves 
"are normally expected. In this paper, we show that, owing to the appearance of traveling waves, the primary Bjerknes 
force near the emitter becomes very large and tends to expel the bubbles up to a stagnation point. Two-dimensional 
axi-symmetric computations of the acoustic field created by a large area immersed sonotrode are also performed, and the 
^^paths of the bubbles in the resulting Bjerknes force field are sketched. Cone bubble structures are recovered and compare 
'reasonably well to reported experimental results. The underlying mechanisms yielding such structures is examined, and 
' it is found that the conical structure is generic and results from the appearance a sound velocity gradient along the 
transducer area. Finally, a more complex system, similar to an ultrasonic bath, in which the sound field results from 
the flexural vibrations of a thin plate, is also simulated. The calculated bubble paths reveal the appearance of other 
^ commonly observed structures in such configurations, such as streamers and flare structures. 

C/3 JKeywords: Acoustic cavitation. Bubble structures. Cavitation fields. Ultrasonic reactors 
>:^-PACS: 43.25.Yw, 43.35.Ei, 43.25. Gf 
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1. Introduction 

A common observation in acoustic cavitation experi- 
ments is the rapid translational motion of the bubbles rel- 
ative to the liquid, and their self-organization into various 
spectacular structures. These structures have been sys- 
tematically reviewed recently and some of them have 
been successfully explained by results derived from single 
bubble physics [l|-|3|. 
I The origin for bubble translational motion in an acous- 
tic field is the so-called Bjerknes force which is the 
average over one oscillation period of the generalized buoy- 
ancy force exerting on any body in an accelerating liq- 
uid . It is commonly expressed in terms of the pressure 
■gradient as: 

Fb = - {VVp) , (1) 

where (.) denotes the average over one acoustic period, 
V is the bubble volume and p the acoustic pressure which 
would exist in the liquid at the center of the bubble if the 
latter were not present. Since V and p are oscillatory, the 
average of their product can be non-zero. 



The most spectacular and known manifestation of the 
primary Bjerknes force occurs in standing waves. It can 
be simply deduced from linear theory that bubbles smaller 
than the resonant size are attracted by pressure antinodes, 
whereas bubbles larger than resonant size are attracted by 
pressure nodes [1] . This was confirmed by the early exper- 
iments of Crum & EUer Q, and is the basic principle of 
levitation experiments used to study single bubble sono- 
luminescence, where attraction by the central antinode of 
the flask counteracts the buoyancy force l^, U | . 



However, nonlinear effects can produce repulsion of 
inertial bubbles from pressure antinodes above a given 
threshold Q . This threshold can be estimated analytically 
for low frequency driving and is found to be near 170 kPa, 
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with a slight dependence on surface tension [12|. This 
can be evidenced in multi-bubble experiments by a void 
region near the pressure antinode surrounded by bubbles 
accumulating near the threshold zone 0, or by bubbles 
self-arrangement into parallel layers shifted relative to the 
antinodal planes [l[, correctly predicted by particle model 
simulations. 

A more important issue concerns the Bjerknes force 
exerted on bubbles by large amplitude traveling waves. 
While small amplitude traveling waves exert a negligible 
Bjerknes force on bubbles, this is no longer true for large 
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traveling waves, and theory predicts a large Bjerknes force, 
oriented in the direction of the wave propagation [J, 0, [l^ . 
This means that an ultrasonic source emitting a travel- 
ing wave would strongly repel the bubbles nucleating on 
its surface. This issue has been investigated theoretically 
by Koch and co-workers [l^ . who, assuming an arbitrary 
wave, traveling in the sonotrode direction and standing in 
the perpendicular plane, showed that the conical bubble 
structure observed under large area transducers [isl . 16 1 
could be partially reproduced by particle models. 

Other complex bubble structures can be observed in 
other configurations, such as ultrasonic baths, and were 
conjectured to result from a combination of both traveling 
and standing waves [l[ . This raises the issue of the origin 
of such traveling waves, which was one of the motivation of 
the present paper and the companion one (which we will 
denote hereinafter by [OL I]). In the latter, we showed 
that traveling waves appear as a simple consequence of 
the attenuation by inertial bubbles. The model presented 
in [OL I] constitutes therefore the missing link in the the- 
ory, and allows to calculate the acoustic field without any 
a priori on its structure, just from the knowledge of the vi- 
brations of the ultrasonic emitter, whatever its complexity. 
From there, the Bjerknes force field can be calculated, and 
the shape of the structures formed by the bubble paths in 
the liquid can be examined. 

Before going further, one should remind that Eq. (IlJ 
is an over-simplification of the complex problem of bub- 
bles translational motion. A correct representation of the 
latter requires to write Newton's second law for the bub- 
ble, accounting not only for the instantaneous driving force 
—Wp [of which ([Ij is the time-average], but also for vis- 
cous drag and added-mass forces . AH members of such 
an equation are dependent of the bubble radial dynamics, 
so that considering Eq. ([T]) as a mean force pushing the 
bubbles is a reduced view of the reality, masking the pe- 
riodic translational motion superimposed to the (macro- 
scopically visible) average translational motion. This is 
historically justified, since the first studies on Bjerknes 
forces aimed at localizing the stagnation points of the bub- 
bles where Fb vanishes 17|, as is the case in the center 
of single-bubble levitation experiments [13, [il. Slightly 
extrapolating this point of view, if one accepts that the av- 
erage bubble velocity can be obtained approximately by a 
balance between Eq. ([T]) and an average viscous drag force, 
a terminal mean velocity of the bubble can be calculated, 
which allowed for example successful particle simulations 
of bubble structures [H . This raises the issue of nontrivial 
averaging procedure for moderate or large drivings 1^ 2^ , 
which may be performed by elaborate multiple scales pro- 



cedures [21|. However, some experimental situations ex- 
ist where such a terminal velocity cannot be defined, and 
a bubble may wander between the nodes and antinodes 
of a standing wave 22|. The description of such a phe- 
nomenon requires the simultaneous resolution of the in- 
stantaneous radial and translational equations of the bub- 
ble, initially proposed in Ref. 23|, and improved recently 



by a Lagrangian formulation 23"2B|- The main result of 
the latter studies is that the radial and translational mo- 
tions are strongly coupled, so that the bubble dynamics 
equation is also affected by the translational motion. Di- 
rect simulation of the two coupled equations in standing 
waves fields reveal that, apart from the classical scheme of 
bubble migration toward stagnation points, some bubbles 
may have no spatial attractors and can wander indefinitely 
between a node and an antinode, as observed in pjj . Such 
a behavior, termed as "translationally unstable" , has been 
found to result from an hysterctic response of the radial 
bubble dynamics below the main resonance (27| . 

In spite of the latter remarks, we will keep in this paper 
the classical picture of the mean primary Bjerknes force 
defined by Eq. ([1} acting on the bubbles, and focus on 
the effects of traveling waves. The paper is organized as 
follows: in section 2, we will first briefly recall the main 
results on the primary Bjerknes force, and indicate how 
it can be calculated for an arbitrary bubble dynamics in 
a given acoustic field. In section 3, we will calculate the 
Bjerknes force field in acoustic fields calculated with the 
model proposed in [OL I], which we will briefly recall in 
section 3.1. First, in section 3.2, the ID configuration ex- 
amined in [OL I] will be considered. Then, in section 3.3, 
we will examine a 2D axi-symmetrical configuration, con- 
stituted by a large area sonotrode emitting in a large bath. 
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similar to the experiments reported in Refs. [15, 
Finally, section 3.4 will address another 2D configuration, 
mimicking an ultrasonic bath in which the acoustic field 
is produced by a plate undergoing fiexural vibrations. For 
both 2D configurations, the bubble paths will be drawn 
from the knowledge of the acoustic and Bjerknes force 
fields at every point in the liquid. The structures obtained 
will be compared to experimental results of the literature 
and discussed. 



2. Primary Bjerknes forces 

2.1. Intuitive analysis and linear case 

The physical origin of the primary Bjerknes force can 
be recalled simply by considering a mono-dimensional wave. 
The instantaneous pressure force exerted by the external 
liquid on a liquid sphere that would replace the bubble 
is approximately the difference Ap between the instan- 
taneous acoustic pressures on two opposite sizes of the 
sphere, multiplied by the bubble area S. Besides, the pres- 
sure difference Ap is roughly dp/dx x 2R, so that the in- 
stantaneous force is roughly 2R x S x dp/dx ~ Vdp/dx. 
Generalizing this result in 3D Eq. ^ is recovered. 

Along an acoustic cycle, the bubble therefore wanders 
forward and backward along the direction of the pressure 
gradient, under the influence of this instantaneous force, 
but the two motions may not exactly compensate, because 
the bubble may be for example larger when the pressure 
gradient is directed forward than when it is directed back- 
ward. The average force is therefore a matter of phase 
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between the volume V and the pressure gradient dp/dx, 
which can be better understood with the schematic repre- 
sentation of Fig. [TJ the phase shift between V and dp/dx 
can be decomposed into two part: the phase shift 7 be- 
tween volume V and pressure p, and the phase shift 9 
between p and its gradient dp/dx. The former depends 
on the way the bubble responds to the local acoustic field, 
that is on the bubble dynamics, whereas the latter depends 
on the structure of the acoustic field. All the results men- 
tioned in the introduction can be interpreted from this 
picture. 



frequency cj, and defined in any point r by 



7 



e 



V 



p 



Bubble 
dynamics 



dp 

dx 



Acoustic field 



Figure 1: Schematic interpretation of Eq. JTJ. The phase shift be- 
tween bubble volume V and pressure gradient dp/dx can be decom- 
posed into two parts: 7 depending on the bubble dynamics, and 9 
depending on the acoustic field. 

Let's take for example the case of sub- resonant bubble 
oscillating linearly. In this case, pressure and volume are 
in opposition (7 = tt). In a pure ID linear standing wave, 
away from pressure nodes or antinodcs, p and dp/dx are 
either in phase (6* = 0), or in phase opposition (0 = tt), de- 
pending on the location relative to the pressure antinode. 
Thus, the phase shift between V and dp/dx is either or 
TT. This yields therefore a large average value for the prod- 
uct {V dp/dx), except in the pressure antinodes and nodes 
where it is zero. Conversely, for a traveling wave, pressure 
and pressure gradient are in quadrature {9 = 7r/2), so that 
(Vdp/dx) is clearly zero in this case. 

If the analysis is rather simple for linear oscillations, 
this is no longer the case for strongly nonlinear inertial os- 
cillations. In that case, the bubble radial motion is mainly 
driven by the inertia of the liquid, and the bubble radius 
contains a large out-of-phase component with respect to 
the driving pressure p. This is the reason why pressure 
antinodes may become repulsive even for a sub-resonant 
bubble in a large amplitude standing wave @ , and why the 
Bjerknes force may become very large in traveling waves 
(isl . . The next section quantifies this qualitative anal- 
ysis 

2.2. General calculation of the Bjerknes force 

The model described in [OL I] was shown to result from 
the assumption that the bubbles mainly respond to the 
first harmonic of the field, which we termed as "first har- 
monic approximation" (FHA). We therefore assume that 
the pressure field in the liquid is mono-harmonic at angular 



p(r,i)=po + i(P(r)e^"* 



P(r)e" 



(2) 



which, writing P = |-P(r)| exp(i0(r)), can be recast as 

p(r,i) =po + |P(r)|cosM + 0(r)] (3) 

This expression may represent a traveling wave, a standing 
wave, or any combination of both. We also define the 
pressure gradient in general form as 



dp 

dxi 



(r,i) =G,(r)cos [wt + V,(r)], 



(4) 



where the fields Gi and ipi can be expressed as functions 
of P and (p once the acoustic field is known. 

The following two extreme cases deserve special consid- 
eration: for a standing wave, (/)(r) = (j>o, so that Gi(r) = 
dP/dxi and V'i(r) = 0o; for a traveling wave, P{y) = Pq 
and ^(r) = — k.r so that Gi(r) = kiP^ and ^i[r) = 

m - ^/2. 

The expression of the Bjerknes force on the bubble lo- 
cated at r reads, from ([T]): 



V{v,t)cos[ujt + iy,{v)] dt, (5) 



where T is the acoustic period and V{r, t) is the instanta- 
neous volume of a bubble located at r and can be calcu- 
lated by solving a radial dynamics equation, for example: 



2ct , P , , 
- - -p(r,t). (6) 



Pi[RR+ -R' 



The bubble volume depends on r because two bubbles lo- 
cated at different points may be excited by fields of differ- 
ent amplitudes |P| but also different phases (j). However, 
in order to be able to calculate the volume V of any bubble 
over one acoustic period independently of its spatial loca- 
tion, we must fix the phase of the driving field in Eq. ^ by 
a convenient change of variables. We therefore set, taking 
this opportunity to non-dimensionalize the variables: 

p(r,t) =po(l- |P*|cosT*), (7) 

where the minus sign has been chosen to be consistent with 



earlier studies |12l . l29j . and, comparing this expression with 
Eq. dSl), we get: 



P* = P(r)/po, 
r*(r,t) =ujt + 4){y) ~ TT. 



(8) 
(9) 



We now note V{t*) the volume of the bubble when it is 
driven by the pressure field ([7]) , and making the change of 
variables in ([S]), we get: 



G.W— ]^ F(r*)cos[r* 



(10) 
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which we recast as: 

Fb,: = G,(r)(/ccos[(/)(r)-i^,(r)] 



with 



+/ssin [0(r) - tlji{r)] 



1 r^^ 

Ic^^ V{t*)cost* dT*, 
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1 f^^ 

— / V (t* ) s.m T* dr* . 
2^io 



(11) 

(12) 
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This decomposition (which was suggested in a shghtly dif- 
ferent form by Mettin []]|), has the advantage to clearly 
decompose the respective influences of the bubble dynam- 
ics through integrals Ic and Is, and the one of the acous- 
tic field, through the phase shift 0(r) — V'i(r)- In the case 
of a ID wave, the latter corresponds to the angle 6 that 
we defined in figure [1] Integral Ic measures the standing 
wave contribution to the Bjcrkncs force (0 — t/" = or tt), 
while integral Is measures the traveling wave contribution 
(0- V = ±7r/2). 

The two integrals Ic and Is can be easily calculated 
numerically by solving a radial dynamics equation to ob- 
tain V and averaging over one period. They can also be 
calculated analytically, trivially for linear oscillations, and 
in a more complex manner for inertial oscillations in a 
small size range above the Blake threshold [l^l- The cor- 
rect matching between the two extreme cases is however 
difficult, so that we will use the numerical values here- 
inafter. 

The results are displayed in Fig. [2] in the case of air 
bubbles in water in ambient conditions for two ambient 
radii i?o = 3 /im (solid lines) and i?o = 5 /im (dashed 
lines) . The integral Ic is represented in signed logarithmic 
scale (note the different scales for the positive and nega- 
tive parts). It is seen that it is positive for low drivings, 
and quickly increases near the Blake threshold by about 3 
orders of magnitude. In this range of acoustic pressures, 
antinodes arc attractive. Then, above \P*\ = 1.7 bar, 1^ 
becomes largely negative and the antinodes become repul- 
sive, as was found in Ref. Q . 

The integral Is is always positive, is very weak below 
the Blake threshold, and drastically increases above the 
Blake threshold (by 6 to 7 orders of magnitude). This 
predicts that, as reported in Ref. [l^l, the Bjerknes force 
can become very large in traveling waves. 

3. Results 

3.1. Simulation method 

The complex acoustic field P is obtained by solving a 
nonlinear Helmholtz equation, which has been detailed in 
[OL I] and is briefiy recalled here for completeness: 
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Figure 2: Average quantities Ic and Ig for an air bubble of ambient 
radius 3 fim (solid line) and 5 fira (dashed line) in water, as a function 
of the dimensionless acoustic pressure |P*| for frequency of 20 kHz. 
The vertical lines represent the Blake thresholds for 3 fim (solid) and 
5 bubbles (dashed). 



where the complex wave number is given by: 
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The bubble number N is defined as a step function: it is 
assumed zero in the zones where the acoustic pressure is 
less than the Blake threshold, and is assigned to a constant 
value in the opposite case. 
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ff \P\ > Pb 
ff \P\ < Pb 



(17) 



\7^P + P{\P\)P ^0. 
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3.2. ID results 

We first calculate the primary Bjerknes force for the 
same ID configuration as in [OL I]: the domain length is 
10 cm, the right boundary is assumed infinitely soft, air 
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bubbles of ambient radius 5 /im in water. For a low ampli- 
tude of the emitter {Uq ~ 0.3 /im), the pressure amplitude 
profile is almost a perfect standing wave (Fig. [5^), and 
Fig. ^jp exhibits the classical picture of a somewhat low 
primary Bjerknes force pushing the bubbles towards the 
antinodes (note that the distortion of the force profile in 
Fig. [SJd is an artifact of the signed logarithmic scale used 
in ordinate). In this case, the value of the Bjerknes force 
mainly owes to the Ic term in Eq. (|lip . 
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Figure 3: Simulation of the sound field and resulting primary Bjerk- 
nes force for 5 /im air bubbles in water, with an emitter displacement 
amplitude is Uo= 0.3 fira at 20 kHz. (a): Pressure amplitude profile. 
The horizontal dashed line represents the Blake threshold for 5 fira 
bubbles, (b): Bjerknes force exerted on bubbles at each point of the 
domain (with signed logarithmic scale in ordinate). The circles in 
Fig. (a) represents the stable stagnation points for the bubble. 

For larger emitter amplitude {Uq = 5 ^tm), the pressure 
amplitude profile is strongly damped near the emitter, as 
already commented in [OL I] (Fig. |4^). This strong at- 
tenuation produces a noticeable traveling part in the wave 
and thus a large 1$ term in Eq. which, as shown in 

Fig. Hb, results in a huge positive force near the emitter 
almost 6 orders of magnitude higher than the maximal 
force visible in Fig |3Jd (note that the scales used in both 
figures are identical). The bubbles would consequently be 
strongly expelled from the emitter, and travel right to the 
first stagnation point which is located somewhat far from 
the emitter (see the leftmost circle marker on Fig.HJi). The 
occurrence of such stagnation points has been proposed to 
be the key mechanism for bubble conical structure forma- 
tion [l|, [iJI . We will show in the next section that this is 
partially true, but involves some additional subtleties. 
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Figure 4: Simulation of the sound field and resulting primary Bjerk- 
nes force for 5 fim air bubbles in water, with an emitter displacement 
amplitude is t/o= 5 /im at 20 kHz. (a): Pressure amplitude profile. 
The horizontal dashed line represents the Blake threshold for 5 fj,m 
bubbles, (b): Bjerknes force exerted on bubbles at each point of the 
domain (with signed logarithmic scale in ordinate). The circles in 
Fig. (a) represents the stable stagnation points for the bubble. 



3. 3. Sonotrode 

3.3.1. Experiments with a 12 cm diameter sonotrode 

In this section we compare the results of our model to 
cone bubble structures images presented in Ref. jl5| . The 
latter work used a 100 cm x 60 cm rectangular water tank 
of 40 cm depth, where a sonotrode of diameter 2a = 12 cm 
driven at 20.7 kHz is immersed at 3 cm below the free 
hquid level (Fig. [5]). 

The geometry simulated follows the experimental con- 
figuration used in Ref. [l3| as closely as possible. However, 
in order to avoid time-consuming 3D simulations, we re- 
placed the rectangular tank by a cylindrical one with the 
same depth and a diameter of 60 cm, and simulate only a 
half-plane cut in axi-symmetrical mode. The characteris- 
tic of the bottom of the tank is not specified in Ref. [isj 
but following a similar work by the same authors in which 
the same bubble structures were observed [lij . we con- 
sidered an anechoic tank bottom. The lateral sides of the 
tank were taken as infinitely rigid boundaries, and the free 
liquid surface as infinitely soft. The nonlinear Helmholtz 
equation with the latter boundary conditions was solved 
in axi-symmetrical geometry, using the commercial COM- 
SOL software. 

The transducer is also simulated in order to account 
for its lateral deformation in the liquid, which, as will be 
seen below, is necessary to catch some experimental fea- 
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Figure 5: Axi-symmctrical geometry for large area sonotrode. 

tures. However, since details on the internal structure of 
the transducer are not given in Refs. [H, [IBl, we assumed 
the latter made of steel and following a non-dissipative 
elastic behavior represented by Hooke's law (with Young 
modulus E = 2 X 10^^ Pa, Poisson ratio f = 0.3, den- 
sity ps = 7900 kg/m'^). The vibration of the transducer 
is coupled to the acoustic field in the liquid by using the 
convenient cinematic and dynamic interface conditions, as 
detailed in Ref. [s^] ■ We simulate only the bottom part of 
the transducer, containing the whole part immersed in the 
liquid and an arbitrary small length (3 cm) of the emerged 
part. A uniform sinusoidal displacement of amplitude Uq 
is imposed on the upper boundary of the simulated trans- 
ducer rod. In order to match the conditions of the ex- 
periments, the acoustic intensity / entering the medium 
through the lower boundary of the sonotrode is calculated 
by: M 

s 

Both parameters Uq and iVo are varied in order to obtain 
the required value of /. 

In order to tentatively exhibit the bubble structures 
formed in a given configuration, the bubble paths, gener- 
ally termed as "streamers" , will be materialized by draw- 
ing the "streamlines" of the Bjerkncs force field in some 
parts of the liquid. The adequate choice of the starting 
points of these streamlines is difficult, because it would 
require a clear knowledge of the bubble nucleation pro- 
cess. Solid boundaries are known to act as sources of 
bubble nuclei, where the latter may be trapped by mi- 
croscopic crevices (sH . Common observation of cavitation 
experiments indeed show that bubbles often originate from 
the transducer area, which might suggest that the release 
of crevice-trapped bubbles is more efficient on vibrating 
surfaces. We will therefore launch systematically stream- 
lines of the Bjerknes force field from equidistant points 
of vibrating boundaries, and we will term them as "S- 



streamers" . 

Besides, many bubble structures appear far from solid 
boundaries as a more or less complex set of bubble fila- 
ments In that case, bubble seem to originate from 
given points of the bulk liquid, but the precise mechanism 
of nucleation of such bubbles is not clear. Although it has 
long been thought that sub-micronic nuclei could grow up 
to the Blake threshold by rectified diffusion [H, [sl] , this 
is ruled out by nonlinear theory, since a sub-Blake bub- 
ble cannot grow by rectified diffusion [s^l- Coalescence is 
an alternative growth process, but this issue is yet unre- 
solved. We will therefore assume that a bubble is visible 
and contributes to structures only if it is inertially oscil- 
lating, that is in zones above the Blake threshold. This 
is anyway consistent with our assumption on the bubble 
density ([TT]) used to calculate the acoustic field. We will 
therefore launch streamlines from arbitrary points located 
on the calculated curves \P\ = Pb, where Pb is the Blake 
threshold {Pb = 1.178 bar for 2 ^m bubbles in ambient 
conditions, sec [OL I]). We will refer to such streamlines 
as "L-streamers" and we will represent then with a color 
different from surface streamlines in order to distinguish 
them. 

Figure El displays one of the original images obtained 
in Ref. [l5| for an acoustic intensity / = 8.2 W. The 
cone is completely formed and ends in a long tail under- 
going lateral fluctuations, which explains the slightly non- 
symmetric shape of the structure. Besides, it can be easily 
seen that a large region inside the cone is poorly populated 
in bubbles, compared to the the immediate vicinity of the 
transducer and the lateral boundaries of the cone. 




Figure 6: Original image of cone bubble structure (reprinted from 
Ref. [H, with permission from Elsevier). The transducer diameter 
is 12 cm. 

We present in Fig. [7] a comparison between this pic- 
ture and the result of our model for Uq = 1.4 fim and 
A'o = 360 bubbles. mm^"^. Because the original picture is 
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non-symmetric, and to ensure a maximal objectivity in the 
comparison, we present our result compared both to the 
left part of the cone (Fig. [7^) and to its mirrored right part 
(Fig. Uh)- Besides, the original picture has been video- 
reversed in order to obtain black bubble paths on a white 
background. We emphasize that this was the only image 
treatment performed. The Blake threshold contour curve 
is displayed in thick solid red line. The S-streamers, orig- 
inating from the transducer, are displayed in black, while 
the L-strcamcrs. originating from arbitrary points on the 
Blake threshold contour line are displayed in blue. 
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Figure 7: Comparison between the experimental picture of Fig. |6] 
(the image of Ref [Tsll has been video reversed in order to make the 
comparison easier) and the numerical results obtained for 2 fim air 
bubbles, with Uq = 1.4 fim, Nq = 360 bubbles. mm"''. The input 
intensity is 8.2 W.cm"'^. The black lines arc the S-streamers, the 
blue lines are the L-streamers, and the thick solid red line is the 
Blake threshold contour curve {\P*\ = 1.178 for 2 fim bubbles), (a) 
comparison with the left part of Fig. [6] ; (b) comparison with the 
mirrored right part of Fig. \E\ The horizontal and vertical lines in the 
liquid just mark the separation between various subdomains and do 
not have any physical meaning. 



First, it is seen that the global shape of the cone is 
correctly reproduced. The acoustic pressure on the axial 
point of the emitter is 2.16 bar. Cavitation occurring near 
the sonotrode dissipates a lot of energy, which produces 
a strong attenuation and therefore a large traveling wave 
contribution in the vertical direction. Bubbles originating 
from the transducer (black lines in Fig. [7]) are therefore 
strongly expelled from the sonotrode surface, the mecha- 
nism being the same as the one explained above for ID 
waves (see Fig. 2]). 

A more detailed analysis of the vertical component of 
the Bjerknes force can be made by referring to the z- 
projection of Eq. ([TT|). The green line in Fig. [5] represents 
the standing wave contribution Ic cos [4>{r) — V'z(r)] as a 
function of the distance to the sonotrode (the latter be- 
ing located on the right of the graph), the red line is the 
traveling wave contribution sin [0(r) — '(/'^(r)], and the 
blue line is the sum of the two latter. The sign of the 
blue line represents therefore the sign of the z-componcnt 
of Fb, which, if negative, corresponds to a downward ori- 
ented force. It can be seen that near the sonotrode, the 
vertical Bjerknes force is dominated by the strong repul- 
sive traveling wave contribution, as was the case for the 
ID simulation, which is clearly due to the strong attenu- 
ation of the wave near the sonotrode. It can be noticed 
by the way that the standing wave contribution is repul- 
sive only in a small layer near the sonotrode because the 
acoustic pressure is larger than the threshold 1.7 bar in 
this zone. Then, slightly before z = —0.01 m, the positive 
standing wave contribution cancels exactly the traveling 
wave one, and becomes dominant, so that the Bjerknes 
force becomes positive. This means that, as far as only 
the z-component of the Bjerknes force is concerned, there 
is a stagnation point for bubbles near z = —0.01 m. In fact 
as will be seen below, this point is a saddle-point since the 
axis is repulsive in the radial direction at this point, so 
that the bubbles expelled from the sonotrode brake here 
in the z-direction and follow their motion radially, which 
explains the formation of the void region in the core of the 
cone. Near z = —0.02 m, the standing wave contribution 
changes sign again, so that the force becomes downwards 
and dominated by the traveling wave again. Drawing the 
same curves for larger distances from the sonotrode would 
show that this is the case up to another sign-change of the 
standing wave contribution, near z = —0.155 m, which 
constitutes a real stagnation point since there, the radial 
force is oriented towards the axis. Contrarily to earlier in- 
terpretations [l|, [3] , the present results suggest that the 
tip of the cone (near z = —0.07 m) is not a stagnation 
point, and that the real one is located well below, so that 
the bubbles follow closely the axis up to the latter on an 
appreciable distance. This is consistent, at least qualita- 
tively, with the original picture of the cone Fig. which 
shows that the cone ends into a long fluctuating tail. 

We now look at the behavior in the radial direction. 
Figure ini displays sin^6'r, where Or = (/)(r) — ^/'r(r) is the 
phase shift between p and dp/dr. It can be seen that a 
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Figure 8: Blue line: magnitude of the parenthesis of Fbz in Eq. {TTJ 
along the symmetry axis. The emitter is on the right of the graph. 
Green line: cos term in the parenthesis of Fb^ in Eq. Illll l (standing 
wave contribution). Red line; sin term in the parenthesis of Fbz in 
Eq. Illlll (traveling wave contribution). 



large region of traveling wave in the r-direction (sin^ 9r — 
1) surrounds approximately the cone boundary, but that 
sin^ 9r decreases back to zero when either entering the core 
of the cone, or moving outward perpendicular to the cone 
boundary. In the latter two regions, the wave has there- 
fore a larger standing part, so that the classical picture 
of attraction by pressure antinodes and repulsion by pres- 
sure nodes applies. Thus, when the pressure is maximal 
on the axis, bubbles converge toward the latter, and this 
explains the formation of the narrow cone tip. This is the 
case for z = —0.04 m (thin solid line in Fig. [TU)) . The op- 
posite holds in the core of the cone, where the variations 
of the acoustic pressure in the radial direction presents a 
local minimum on the axis, for example at z = —0.01 m 
(dashed line in Fig.ITU)) or z = —0.02 m (dash-dotted line). 
The radial component of the Bjerknes force in this zone 
is therefore oriented outwards. This is why, as mentioned 
above, the point on the axis (near z = —0.01 m) where the 
z-component of the Bjerknes force change sign is in fact a 
saddle-point, which locally pushes the bubbles far from the 
axis, and produces a void region in the heart of the cone, 
clearly visible on the experimental picture. This feature 
has been commented in Ref. [isj and was attributed to the 
nonlinear reversal of the Bjerknes force in standing waves 
near 1.7 bar. Our results suggest that this is not the case, 
and that the void region results from a combination of a 
canceling z-component of the Bjerknes force and a local 
inversion of the radial standing wave pressure profile. 

The shape of our predicted void region shows reason- 
able agreement with the experiments. Furthermore, the 
experimental cone tip seems to be more dense than its core. 
This may be due to nucleation of bubbles at the Blake 
threshold in this part, as suggested by the L-streamers 
starting from the Blake contour loop just above the cone 
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Figure 9: Color plot of sin^ 9r, where 9r represents the phase between 
p and dp /dr. The S-streamers are recalled in black lines. 



tip in Fig. [7] (blue lines). 

Another common observation on such sonotrodes is the 
presence of small streamers on their lateral side, visible 
near the upper left corner of Fig. [7}d. As seen in the right 
part of the latter figure, this phenomenon is reasonably 
caught by the simulation, and this is the reason why the 
deformation of the transducer was accounted for. Indeed, 
our result suggests that such small structures result from 
the lateral vibration of the sonotrode, which emits a radial 
wave, and produces a small zone of large acoustic pressure. 
The bubbles in this zone strongly attenuate the wave, and 
produces a traveling part in the radial wave. The physical 
mechanism is therefore similar to the cone formation, but 
here the stagnation point is very close to the sonotrode 
surface, so that only a small flat filamentary structure is 
formed. 

Figure [TT] shows the same result as Fig. [7}d where we 
sketched additional streamers originating from arbitrary 
points in the liquid (green lines), which makes the com- 
parison of the lateral filamentary structure with experi- 
ments more striking. Furthermore, this representation al- 
lows to evidence streamers starting near the cone lateral 
boundary and quickly merging with the latter, as indeed 
visible on the experimental picture. Figure [TT] also shows 
that the corner of the sonotrode acts as a separatrix be- 
tween the streamers attracted by the cone and the ones at- 
tracted by the lateral filamentary structure, as can be also 
speculated from the experimental picture. We note how- 
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Figure 10: Radial pressure profiles at 2 = m (thick solid line), z = 
—0.01 m(dashed line), z = -0.02 m (dash-dotted line), z = —0.04 m 
(thin solid line). 



ever that streamers starting at a larger distance from the 
cone (see magenta hnes in Fig. [11]) start upwards, whereas 
such streamers on the experimental picture seems to start 
downwards. The experimental image also suggests that 
some streamers starting from points far from the cone (say, 
near r = 0.06 m, z = —0.07 m) seem to be attracted by 
points located outside the picture, and this feature is not 
caught by our simulation. 
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Figure 11: Same as Fig. [TJ^, sketching additional streamers originat- 
ing from arbitrary points points below the Blake threshold (green 
and magenta lines). 

The interpretation of the cone structure as the result 
of the combination of a longitudinal traveling wave and a 
lateral standing wave proposed in Ref. [3| is therefore con- 
firmed by the present model. However, the above analysis 
does not tell much about how such an acoustic field ap- 
pears. Since cone bubble structure arc very robust against 
amplitude, sonotrodc size (see discussion in Ref. [15|), and 



even appear near the walls of ultrasonic baths (see Ref. [ij , 
and next section), there must be some generic mechanism 
responsible of its formation. Since traveling waves consti- 



tute the key phenomenon of the problem, it is interesting 
to sketch the contour lines of the acoustic field phase (0 in 
Eq. ([3])). For a traveling wave these lines are orthogonal 
to the direction of propagation and constitute therefore a 
powerful visual tool to assess the latter. The result is dis- 
played in Fig. [12] (white lines). It can be seen that, while 
the wave mainly propagates along the z-direction inside 
the cone, it bends into an oblique direction near the cone 
boundary, targeting at a point located on the symmetry 
axis (S-streamcrs arc recalled in black). We emphasize 
that the emission of an oblique wave from the sonotrode 
corner was found to occur in all our simulations, whatever 
the sizes of the sonotrode and the liquid domain, and the 
type of the bottom liquid boundary. More importantly, 
we found that the same phenomenon occurs whenever the 
deformation of the sonotrode was accounted for or not, 
which rules out any effect of a non uniform displacement 
of the sonotrode tip. 
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Figure 12: Color plot of the acoustic pressure field in the liquid. 
The white lines are the contour lines of the phase <p of the acoustic 
pressure field (see Eq. The S-streamers are recalled by the black 

lines. 

We therefore infer that the slanting of the wave prop- 
agation direction is definitely linked to the presence of 
strongly driven bubbles near the vibrating area. These 
bubbles dissipate a lot of energy, rendering the square 
of the local wave number almost purely imaginary. But 
the acoustic field near the outer points of the sonotrode is 
weaker, because this region is less constrained laterally (see 
pressure profile in thick solid line in Fig. fTU)) . Therefore, as 
evidenced in [OL I] (see Fig. 4 in the latter reference), the 
real part of k is higher in the central part of the sonotrode 
than in its outer part, and the opposite holds for the sound 
velocity uj/^{k), which is confirmed by Fig. [T31 There ap- 
pears therefore an outward gradient of sound speed along 
the sonotrode area, which, following Huygens principle, 
bends the wave number towards the axis, and produces 
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the conical traveling wave visible on Fig. [121 This travel- 
ing wave produces in turn a strong Bjerknes force directed 
along the propagation direction of the wave (a large Is 
term in Eq. (jlip ). which structures the bubbles into a con- 
ical shape. The latter scenario was qualitatively checked 
by simple linear acoustics simulations: setting the sound 
field to ci uniformly in the liquid except in a thin cylinder- 
shaped region below the transducer, the bending of the 
iso-d> lines was indeed observed. 
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Figure 13: Radial profile of sound velocity on the transducer area. 
The sound velocity is non-dimensionalized by the sound velocity in 
the pure liquid c;. 

A remarkable feature of cone bubble structures is the 
invariance of their shape when intensity is increased above 
a given level [l^. Following the suggestion of an anony- 
mous reviewer, we performed an additional simulation of 
the above configuration, increasing the sonotrode displace- 
ment Uq up to 4.2 fim (instead of 1.4 /im). The results are 
presented as supplementary material. The cone shape ob- 
tained is indiscernible from the one presented in Fig. 1111 
However, a close examination of the axial pressure pro- 
files in the two cases reveals that the acoustic fields differ 
mainly in a thin layer of about 5 mm near the sonotrode 
(see supplementary material). This is a clear manifesta- 
tion of the self-saturation effect inherent to the present 
model through the field-dependence of the attenuation co- 
efficient (see [OL I]): increasing the sonotrode displace- 
ment produces a large increase of acoustic pressure only 
locally, but the bubbles in this zone being excited more 
strongly, they dissipate the excess acoustic energy very 
rapidly. It should be noted that experimental manifesta- 
tions of this phenomenon have been reported in the early 
work of Rozenberg [s^ . 

Other studies demonstrates that conversely, for low ex- 
citations, the shape of the cone bubble structure does de- 
pend on the driving. Although some of our simulations 
could partially catch such a dependence, convergence prob- 
lems in this range prohibited any firm conclusion. We ob- 
served however that the cone shape was very sensitive to 
the choice of the bubble density when the latter was low 
enough. This might suggest that the experimentally ob- 
served shape dependence for low drivings would be rather 
due to the variation of the bubble density with acoustic 
pressure, than to the pressure dependence of a single bub- 



ble dissipation. This is a missing brick in our model since 
we consider constant bubble densities above the Blake 
threshold. We thus infer that the model in its present 
form is not able to catch the latter experimental feature. 

One should finally mention, as underlined in Ref. [28j |. 
another explanation for the cone structure robustness against 
driving level above a certain threshold, borrowed to phase 
transitions theory. Skokov and co-workers measured laser 
intensity transmission through a cavitation zone and ob- 
tained time-series presenting fiuctuations whose power spec- 
trum were found to be inversely proportional to frequency 
[36l . I37I I . This feature, also termed as "flicker noise" has 
been shown to occur in various physical processes (ssj . 
and has been interpreted as a consequence of the interac- 
tion between two phase transitions, one subcritical, the 
other supercritical [s^. The result is that the system 
self-organizes into a critical state, whatever the precise 
value of the controlling parameter, contrarily to classi- 
cal critical states which require a fine tuning of the lat- 
ter to be reached. This phenomenon has been termed as 
"self-criticality" and produces organized self-similar spa- 
tial structures, reminiscent of bubble web-like organiza- 
tion. Self-criticality can be modeled generically by a stochas- 
tic dynamical system of two equations, generalizing the 
Ginzburg-Landau equation. One of the two order param- 
eters exhibits fluctuations between two attractors, with 
a power spectrum varying in 1//, as observed in experi- 
ments. This original theory has the advantage to explain 
some generally overlooked features of cavitation clouds by 
a universal physical mechanism. However it still remains 
very far from the precise cavitation physics, and the pro- 
posed equations are phenomenological. In particular the 
precise physical sense of the order parameters remains to 
be explicited, maybe on the light of the coupled evolutions 
of the bubble field and the acoustic wave. For now, the 
model is still far from a predictive tool for acoustic cavi- 
tation, but this promising approach remains opened. 

3.3.2. Experiments with a 7 cm diameter sonotrode 



Reference jl5| also present results for thinner sonotrodes, 
but the latters produce acoustic currents which deform 
the cone structure, so that a direct comparison of the ex- 
perimental images and simulated cone structures is not 
possible in this case. In spite of the latter restriction, 
we performed additional simulations for a 8 mm diame- 
ter sonotrode (referred as "type B" in Ref. [l3|): h^ an 
otherwise identical geometry, for 2 /im bubble radii and a 
bubble density A'o = 90 bubbles/mm'^, in order to match 
at best the experimental shape of the cone structure. The 
result of the visual comparison between the simulated cone 
structure and Fig. 2 in [l5| is deferred to supplementary 
material, and shows that, in spite of the blurring of the 
structure by acoustic currents, some similarities in the 
cone shape can be observed. 

Apart from imaging cone structures, Dubus and co- 
workers have collected valuable quantitative experimen- 
tal informations in additional studies, using 7 cm diam- 
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Figure 14: Pressure profile on the symmetry axis for a 7 cm diameter 
sonotrode, with Rq = 2 fim, No = 90 bubbles/mm'^. The velocity of 
the sonotrode tip is 1.31 m/s. Solid line: prediction by the present 
model. Stars: experimental measurements redrawn from Fig. 6 in 
Ref. [H. 



eter sonotrodes excited in pulsed mode in order to avoid 
acoustic currents [3 22 1. In order to put our model to 



the test, we performed additional simulations of a 7 cm 
sonotrode, again with for 2 /im bubble radii and A^o = 90 
bubbles/mm'^. The transducer displacement was set in or- 
der to match the velocity of the sonotrode area measured 
in Ref. [l^ (1-31 m/s). It is interesting to compare the 
calculated and experimental axial acoustic pressure pro- 
files (Fig 6 in Ref. [l^, star signs). The result is displayed 
in Fig. [TH It can be seen that the agreement is somewhat 
poor: even if the order of magnitude of the predicted pres- 
sure field is reasonable, it vanishes much more slowly than 
the experimental one. We increased the bubble density up 
to Nq = 360 bubbles/mm^ without appreciable change in 
the pressure profile, except near the sonotrode. This un- 
derlines the limit of the present model, and may result of 
the rough assumption of a constant bubble density. Be- 
sides, the pressure field calculated in the immediate vicin- 
ity of the transducer {z = 0) is found to be much larger 
than the experimental one. It should however be noted 
that hydrophones have a finite size, which does not allow 
to measure the acoustic pressure exactly near 2 = 0, where 
the present theory predicts a large pressure gradient. 

Finally, Dubus and co-workers [28| proposed an alter- 
native explanation of the bubble cone structure. Their 
analysis relies on the formation of a nonlinear resonant 
layer of bubbly liquid attached on the transducer, the fo- 
cusing being qualitatively attributed to the curvature of 
the bubble layer. This layer produces a phase shift of 
the wave emitted by the sonotrode, dependent on its local 
width. The argumentation is supported by phase mea- 
surements of the pressure field in a horizontal plane below 
the layer. We have therefore reported the corresponding 
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Figure 15: Radial variation of the phase of the pressure field, 
at 1.5 cm from the transducer. The quantity represented is 
{<p{0, 1.5 cm) — (f>{r, 1.5 cm)) / {2n), where (f>{r, z) = arg(p(r, z)). The 
solid line represents the theoretical prediction, and the stars are ex- 
perimental results redrawn from Fig. 10 in Ref. [Ssll . The simulations 
parameter are the same as for Fig. 1141 



experimental values, along with the prediction of our sim- 
ulation in Fig. [15] (the simulation parameters are the same 
as for Fig. [T4| . It can be seen that the agreement is very 
good, and the result was found almost unsensitive to an in- 
crease of the bubble density up to iVg = 360 bubbles / mm^ . 

The present calculation shows therefore that the exis- 
tence of a curved resonant layer of bubble is not neces- 
sary to explain the experimental results, although nonlin- 
ear resonant effects cannot be discarded. The wave focuses 
anyway just because the outermost bubbles in the bubble 
layer are more smoothly driven than the central ones, and 
therefore have a lower influence on the local sound speed. 
This a purely nonlinear effect, and is correctly summarized 
by Fig. 4 in [OL I]. Besides, our simulations indeed exhibit 
a bubbly layer near the transducer, but rather than being 
resonant, this layer is in fact found to be very dissipative. 
Such a large dissipation cannot be predicted by using a 
reduced bubble dynamic equation, and is strongly corre- 
lated to the inertial character of the bubble oscillations, as 
shown in our companion paper [OL I] . 

3.4- Cleaning bath 

The system considered is axi-symmetrical and is rep- 
resented in a half-plane cut on Fig. [121 The bottom of 
the bath is a thin circular steel plate of 4 mm thickness 
and 20 cm diameter, which upper side is in contact with 
the liquid, while its lower side is free, except its central 
part (red line on Fig. 1161 r £ [0,3 cm]) which is assumed 
to vibrate with a uniform amplitude Uq = 2 /im and a 
frequency of 20 kHz. This boundary condition is a simpli- 
fication of a system where a piezo-ceramic ring would be 
clamped to the bottom side of the plate and impose an os- 
cillatory displacement. The liquid fills the space above the 
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plate, and is limited laterally by a cylindrical boundary, 
which is assumed infinitely rigid, while the free surface is 
considered as infinitely soft. The plate is assumed elastic 
and its deformation follows the Hooke's law, its vibrations 
being coupled to the acoustic field in the liquid by ade- 
quate interface conditions j30| . 
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Figure 16: Geometry and meshing of the cleaning bath configuration. 
The geometry has axial symmetry and only a half-plane cut is repre- 
sented. The blue line is the symmetry axis. A uniform displacement 
C/q is imposed on the red line. 

Fig. [T7] displays the result obtained with 5 \xa\ air 
bubbles and a water height of 0.14 m. The green line 
represents the deformed shape of the bottom plate, and 
shows that a flexural standing wave is excited in the lat- 
ter. This flexural wave produces a spatially inhomoge- 
neous acoustic pressure field on the plate area, ranging in 
this case between 0.5 and 2.8 bar (the red lines represent 
the Blake threshold contour curves). In the zones on the 
plate where the acoustic pressure is larger than the Blake 
threshold, bubbles are produced, dissipating locally a large 
energy and modifying the sound speed. The remaining 
mechanisms are similar to the cone bubble structure de- 
scribed above. Dissipation produces traveling waves, and 
and can even result into the production of a small cone 
bubble structure (near r = 0.04 m) attached on the plate, 
while other structures on the plate are more similar to the 
above-mentioned streamers attached on the lateral side of 
sonotrodes. A few streamers are also visible in the middle 
of the liquid, located on the antinodes of the wave, which 




Figure 17: Simulation of the cleaning bath configuration. The thick 
red lines materialize the Blake-threshold. The black lines are the S- 
streamers, paths of bubbles originating from the vibrating plate, and 
the blue lines are the L-streamers, path of bubbles originating from 
the Blake threshold. The green thick line represents the deformed 
shape of the plate at ujt = 2n7r. 

An interesting feature can also be seen in Fig. [T7] and 
is magnified in Fig. 118b . It is seen that the S-streamcrs 
(black lines) merging at the vertex of the cone then follow 
a unique path up to a stagnation point located on the sym- 
metry axis. If we also launch bubbles from the Bjerknes 
contour curves, the corresponding L-streamers (blue lines) 
join the main S-streamers up to the stagnation point where 
numerous secondary streamers end up, forming a star-like 
structure. The experimental occurrence of this behavior 
has been reported in a detailed presentation by Mettin [l[ , 
who termed this bubbles arrangement as "flare structure" , 
and is represented in Fig. [18^). It is seen that, apart from 
the curvature of the structure observed in our simulation, 
the latter reproduces reasonably well the main features of 
the phenomenon, especially the cone structure attached 
to the vibrating plate merging into a jet, fed laterally by 
bubbles originating in the bulk liquid. Mettin reported 
that this structure was universally found in cleaning bath 
setups and proposed a qualitative explanation involving a 
"complicated near field structure with shares of both trav- 
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cling and standing waves" . The present results seem to 
enforce this interpretation: a traveling wave takes birth 
near the vibrating area, because of the strong attenuation 
by the bubbles located there, and launches the bubbles far 
from the plate. It ends up at a pressure antinode which 
attracts all the bubbles, either coming from the plate, or 
taking birth in neighboring zones. The lateral enrichment 
of the main bubble path originates from a standing wave 
in the direction perpendicular to the traveling wave. 




1 cm 



Figure 18: Comparison between; (a) a flare structure observed on the 
lateral side of a cleaning bath (reprinted from Ref. [J| with permission 
of Research Signpost; the original image has been 90° rotated for 
comparison purposes) (b) the structure predicted in the simulated 
configuration (zoom of Fig. 1171 with more L-streamers sketched). 



4. Summary and discussion 

The model proposed in the companion paper [OL I] has 
been applied to classical 2D configurations, namely large 
area transducer emitting in a liquid, and cleaning baths. 
The density of bubbles was assumed constant in zones 
where the acoustic pressure is above the Blake threshold, 
and null everywhere else. The possible paths of the bub- 
bles were assumed to originate either from the vibrating 
parts of the solid, or from the Blake threshold contour 
curves, and calculated by computing the primary Bjerk- 
ncs force field directly from nonlinear bubble dynamics 
simulations. 

In the case of large area transducer, cone bubble struc- 
tures observed experimentally can be easily reproduced, 
with reasonable agreement in the cone shape. A strongly 
dissipative bubble layer was found to appear near the trans- 
ducer, and the cone boundaries were shown to be made of 
bubbles following a focused traveling wave. The focusing 
was found to result from a radial acoustic pressure gra- 
dient on the transducer area, which, following the result 
of the companion paper [OL I] produces a radial gradi- 
ent of sound velocity. The streamers located on the lateral 
boundary of the transducer, observed experimentally, were 



correctly reproduced by our simulations, provided that the 
elastic deformations of the emitting transducer were ac- 
counted for. Besides, the calculated stagnation point on 
the symmetry axis was found to be much farther than the 
cone tip, which would explain the long bubble tail visible 
in experimental picture. 

In the cleaning bath configurations, the fiexural vibra- 
tions of the bottom plate were found to produce several 
zones of large acoustic pressures, located near the plate dis- 
placement antinodcs. These high acoustic pressures pro- 
duce locally a thin layer of strongly oscillating bubbles 
dissipating a lot of acoustic energy, which yields damped 
traveling waves and cone-like structures. In some cases, 
the bubbles reaching the cone tip carry on their motion 
along a unique line, ending into a distant pressure antin- 
ode, and laterally enriched by bubbles originating from 
liquid zones excited above the Blake threshold. The ob- 
tained structure is reminiscent of a flare-like structure de- 
scribed in the literature, and known to occur frequently in 
cleaning baths configurations. 

The reasonable success of our model in predicting ab 
initio such structures is encouraging, and seem to show 
that strong energy dissipation by incrtial bubbles is a key 
mechanism ruling the structure of the acoustic field in a 
cavitating medium. It is interesting to note that the self- 
action of the acoustic field evidenced in the present pa- 
per differs from the mechanisms presented in the work of 
Kobelev & Ostrovski In the latter work, the self- 

action of the acoustic field is mediated by its slow influ- 
ence on the bubble population, while here, the mechanism 
is only due to the bubbles radial motion, even for constant 
bubble density. The latter assumption constitutes how- 
ever a weak point of our model, and requires the arbitrary 
choice of two free parameters: the ambient radius of the 
bubbles Rq, and the bubble density, both being assumed 
spatially homogeneous in regions above the Blake thresh- 
old. A more realistic model would require at least the 
spatial redistribution of the bubbles. This may be done 
for example by coupling the nonlinear Helmholtz equa- 
tion used in this paper with a convection-like equation for 
the bubble number density, as done in the linear case in 
Ref. (ioj . which requires the correct estimation of the av- 
erage translational velocity of incrtial bubbles. As already 
mentioned in introduction, this translational motion is de- 
scribed by a somewhat elaborate physics 2^ 27], and the 
estimation of an average velocity raises the complicated is- 
sue of properly averaging the translation equation 1^ 21 1 . 
Along the same line of investigation, an extension of the 
present work could consist in launching some bubbles in 
the acoustic fields presented here, and to calculate their 
paths within the latter by integrating in time the coupled 
equations of radial and translational motion described in 
Refs. [25, 27 1 . Apart from testing the validity of the cur- 
rent approximation, this would also provide a clear picture 
of the dynamics of bubble drift in the studied structures. 
This may reveal some unexpected features, such as bubble 
precession around some points in the liquid, as evidenced 
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in simple standing waves in Ref. |27| . Such refinements 
may be the subject of future work, keeping however as a 
main objective a model simple enough to be used in real 
engineering applications. 

Finally, many other bubble structures can be observed 
experimentally Among the latter, the grouping of bub- 
bles into so-called "clusters" classically appear in numer- 
ous experimental configuration either in the bulk liquid, 
or as hemi-spherical structures near solid boundaries, es- 
pecially in the case of focused ultrasound 41-4^. In some 
aspects, they behave as a single large bubble and may 
collapse as a whole, emitting a complex set of primary 
and secondary shock-waves (see Ref. [i^ for high-speed 
photographs), and yield strongly erosive effects Q- A 
complete theoretical description of such structures is not 
yet available, the difficulty lying in the strong interaction 
between the bubbles constituting the cloud. Clearly, this 
interaction cannot be accounted for by our simple model 
which includes only primary Bjerknes forces. Moreover, 
these structures often present a transitory behavior, mov- 
ing as a whole entity in the liquid, appearing and disap- 
pearing near other structures |4J] , and have even been 
observed as early precursors of conical structures [IHj . As a 
conclusion, we feel therefore that our model in the present 
form cannot account for such structures. 
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t of Figures 

1 Schematic interpretation of Eq. ((!]). The 
phase shift between bubble volume V and 
pressure gradient dp/dx can be decomposed 
into two parts: 7 depending on the bubble 
dynamics, and 9 depending on the acoustic 

field i 

2 Average quantities Ic and Is for an air bub- 
ble of ambient radius 3 fim (solid line) and 
5 /xm (dashed line) in water, as a function of 
the dimensionless acoustic pressure \P*\ for 
frequency of 20 kHz. The vertical lines rep- 
resent the Blake thresholds for 3 /im (solid) 

and 5 ^m bubbles (dashed) 

3 Simulation of the sound field and resulting 
primary Bjerknes force for 5 ^m air bub- 
bles in water, with an emitter displacement 
amplitude is Uq= 0.3 ^m at 20 kHz. (a): 
Pressure amplitude profile. The horizontal 
dashed line represents the Blake threshold 
for 5 /im bubbles, (b): Bjerknes force ex- 
erted on bubbles at each point of the do- 
main (with signed logarithmic scale in ordi- 
nate) . The circles in Fig. (a) represents the 
stable stagnation points for the bubble. . . [H 

4 Simulation of the sound field and resulting 
primary Bjerknes force for 5 ^m air bub- 
bles in water, with an emitter displacement 
amplitude is Uo= 5 fj,m at 20 kHz. (a): 
Pressure amplitude profile. The horizontal 
dashed line represents the Blake threshold 
for 5 /im bubbles, (b): Bjerknes force ex- 
erted on bubbles at each point of the do- 
main (with signed logarithmic scale in ordi- 
nate) . The circles in Fig. (a) represents the 
stable stagnation points for the bubble. . . [HI 

5 Axi-symmctrical geometry for large area sonotrode. 

6 Original image of cone bubble structure (reprinted 
from Ref. [l5|, with permission from Else- 
vier). The transducer diameter is 12 cm. . . |^ 

7 Comparison between the experimental pic- 
ture of Fig. [5] (the image of Ref [l^ has been 
video reversed in order to make the compar- 
ison easier) and the numerical results ob- 
tained for 2 /im air bubbles, with Uq = 
1.4 //m, iVo = 360 bubbles.mm-^. The in- 
put intensity is 8.2 W.cm~^. The black lines 
are the S-streamers, the blue lines are the L- 
strcamers. and the thick solid red line is the 
Blake threshold contour curve {\P*\ = 1.178 
for 2 /im bubbles), (a) comparison with the 
left part of Fig. |6] ; (b) comparison with the 
mirrored right part of Fig. [B] The horizon- 
tal and vertical lines in the liquid just mark 
the separation between various subdomains 

and do not have any physical meaning. ... 
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Blue line: magnitude of the parenthesis of 
in Eq. ([TT|) along the symmetry axis. 
The emitter is on the right of the graph. 
Green line: cos term in the parenthesis of 
Fbz in Eq. ([TT|) (standing wave contribu- 
tion). Red line: sin term in the parenthesis 
of Fbz in Eq. ([TT|) (traveling wave contri- 
bution) i 

Color plot of sin^^r, where Or represents 
the phase between p and dp/dr. The S- 

streamers are recalled in black lines 

Radial pressure profiles at z = m (thick 
solid line), z = —0.01 m(dashed line), z ~ 
-0.02 m (dash-dotted line), z = -0.04 m 

(thin sohd hue) H 

Same as Fig.[7lD, sketching additional stream- 
ers originating from arbitrary points points 
below the Blake threshold (green and ma- 
genta lines) 

Color plot of the acoustic pressure field in 
the liquid. The white lines are the contour 
lines of the phase of the acoustic pressure 
field (see Eq. ([3|)). The S-streamers are re- 
called by the black lines |^ 

Radial profile of sound velocity on the trans- 
ducer area. The sound velocity is non-dimensionalized 
by the sound velocity in the pure liquid c; . 
Pressure profile on the symmetry axis for a 
7 cm diameter sonotrode, with Rq ~ 2 /im, 
iVo = 90 bubbles/mm'^. The velocity of the 
sonotrode tip is 1.31 m/s. Solid line: pre- 
diction by the present model. Stars: exper- 
imental measurements redrawn from Fig. 6 

in Ref. ^ 

Radial variation of the phase of the pres- 
sure field, at 1.5 cm from the transducer. 
The quantity represented is (0(0, 1.5 cm) — 
(/)(r, 1.5 cm))/(27r), where </)(r, z) = arg(p(r, z)) 
The solid line represents the theoretical pre- 
diction, and the stars are experimental re- 
sults redrawn from Fig. 10 in Ref. [2^. The 
simulations parameter are the same as for 



17 Simulation of the cleaning bath configura- 
tion. The thick red lines materialize the 
Blake-threshold. The black lines are the S- 
streamers, paths of bubbles originating from 
the vibrating plate, and the blue lines are 
the L-streamers, path of bubbles originat- 
ing from the Blake threshold. The green 
thick line represents the deformed shape of 
the plate at ujt = 2n-K 

18 Comparison between: (a) a flare structure 
observed on the lateral side of a cleaning 
bath (reprinted from Ref. [3| with permis- 
sion of Research Signpost; the original im- 
age has been 90° rotated for comparison 
purposes) (b) the structure predicted in the 
simulated configuration (zoom of Fig. [T7] 
with more L-streamers sketched) 
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Fig. [H Ill 

Geometry and meshing of the cleaning bath 
configuration. The geometry has axial sym- 
metry and only a half-plane cut is repre- 
sented. The blue line is the symmetry axis. 
A uniform displacement Uq is imposed on 



the red line 12 
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